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ABSTRACT 

We present Spiteer/Infrared Spectrograph (IRS) 5-21 pm spectroscopic maps towards 
12 regions in the Andromeda galaxy (M31). These regions include the nucleus, bulge, 
an active region in the star-forming ring, and 9 other regions chosen to cover a range of 
mid-to-far-infrared colours. In line with previous results, PAH feature ratios (6.2 pm 
and 7.7 pm features compared to the 11.2 pm feature) measured from our extracted 
M31 spectra, except the nucleus, strongly correlate. The equivalent widths of the main 
PAH features, as a function of metallicity and radiation hardness, are consistent with 
those observed for other nearby spiral and starburst galaxies. Reprocessed data from 
the ISOCAM instrument on the Infrared Space Observatory agree with the IRS data; 
early reports of suppressed 6-8 pm features and enhanced 11.3 pm feature intensity 
and FWHM apparently resulted from background-subtraction problems. The nucleus 
does not show any PAH emission but does show strong silicate emission at 9.7 pm. 
Furthermore, different spectral features (11.3 pm PAH emission, silicate emission and 
[Neill] 15.5 pm line emission) have distinct spatial distributions in the nuclear region: 
the silicate emission is strongest towards the stellar nucleus, while the PAH emission 
peaks 15" north of the nucleus. The PAH feature ratios at this position are atypical 
with strong emission at 11.2 pm and 15-20 pm but weak emission at 6-8 pm. The 
nucleus itself is dominated by stellar light giving rise to a strong blue continuum and 
silicate emission. 

Key words: galaxies: individual: M31 - galaxies: ISM - galaxies: nuclei - infrared: 
ISM - ISM: molecules - ISM: lines and bands 


1 INTRODUCTION 


Mid-infrared spectra provide a unique diagnostic tool to un¬ 
derstand the physical conditions in the interstellar medium 
of galaxies. The rich range of spectral features (includ¬ 
ing Polycyclic Aromatic Hydrocarbons (PAHs), atomic fine 
structure lines (e.g. Ne, S) and the amorphous silicate fea¬ 
ture centred at 9.7 pm provide information on dust prop¬ 
erties, radiation field and star formation. With the advent 
of infrared spac e telescopes, such as the Infrared Space Ob¬ 
servatory ( ISO. iKessler et al.l 199d ) and the Spitzer Space 
Telescope (IWerner et al.l 2004 ). mid-infrared spectroscopy 
has become an important method of investigating the in¬ 
frared emission from galaxies. 

PAHs are known as the c arrier of the ubiquitous 
mid-infrared emissio n bands (e.g. lAllamandola et al.lll989l . 
IPuget fc LegerlflOSgl ). They are large hydrocarbon molecules 
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consisting of -^SO-lOO carbon atoms (iTielend l2008ll . The 
main PAH fe atures are seen at 3.3, 6.2, 7.7, 8.6, 11.2 and 
12.7 / rm (e.g.. lGillett et al]ll973l . fGeballe et'^ll985l . |Peetersl 
I2OIII) . These bands are attributed to the vibrational de¬ 
excitation of PAH molecules through b ending and stretch¬ 
ing modes of C-H and C-C bonds (e.g. lAllamandola et aP 
1 19891 . IPuget fc Legerj 1 19891 ). The 6 to 8 micron features 
are thought to originate mostly from ionized PAHs and 
the 3.3 and 11.2 gm emission ban ds from neutral PAHs 
(e.g. Hudgins fc Allamandolal 200^ and reference therein, 
iHonv et alJ 2001 , Galliano et al. 20081 ) . 


The relative strength of the PAH features vary spatially 
within extended ob jects and from galaxy-to-galaxy (e.g. 
iGalliano et al.|[2008l) . While extinction does in fluence the in¬ 
dividual PAH bands to different degrees (e.g. iBrandl et al.l 
I2OO6I : IStock et al.l |201^ . the observed spread in PAH in- 
tensi ty ratios is dominat ed by the PAH charge balance 
(e.g. IGalliano et al.l[2008l ). In addition, PAH intensity ra¬ 
tios change more drastically close to active galactic nuclei 
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where the overall strength of PAH emi s sion also gets weaker 
dRoche et al.lll99i] . ISmith et ahlbOOT^^ . ISmith et alj (l2007bl l 
found that the mid-infrared spectra from some weak AGNs 
show suppressed 6 to 8 fim PAH features (by up to a factor 
of 10 in strength) but are bright at 11.2 /rm. One possible 
explanation for this behaviour is that AGNs alter the grain 
composition by selective destruction of small PAHs and/or 
excite the PAHs. Alternatively, the PAH emission is mod¬ 
ified by the l ow star formation intensity in the centres of 
many AGNs dSmith et al.ll200'^ l. 

Previous studies of nearby galaxies indicate that 
metallicity and radia tion har d ness c an both affe c t the 
PAH emission (e.g. Mj^ldgn 2000|, Beirar)et_alJ l2006l . 
lEngelbracht et aP l2008l . Munoz-Mateos et al.l 2009ll. In a 
sampl e of nearby starburst galaxies lEngelbracht et al.l 
ll2008l f found the PAH equivalent widths (EQWs) to 
decrease with inc reasing radiation hardness; however, 
ISrandl et al.l (|2006l f found no such correlation within their 
starburst sample. In addition, PAH EQWs show an anti¬ 
correlation with metallicity in star-forming galaxies. This 
variation of PAHs among g alaxies has also bee n observed 
within H ii regions of MlOl (iGordon et ahlfioo^ l. But there 
are no other investigations done on a single star-forming 
galaxy with sufficiently high resolution to see whether the 
correlations mentioned above hold within a galaxy similar 
to the Milky Way. 

The amorphous silicate feature at 9.7 fim is another 
aspect of the mid-i nfrared spectra of g alaxies and in par¬ 
ticular their nuclei. ISpoon et al.l ll2007l i classified infrared 
galaxies based on the equivalent width of the 6.2 fim PAH 
feature and the strength of the 9.7 /rm silicate emission or 
absorption feature. They found galaxies spread along two 
distinct branches: one in which silicate absorption strength 
was anti-correlated with PAH equivalent width, and an¬ 
other in which the weak silicate feature strength did not 
depend on the 6.2 /rm equivalent width. Silicate emission 
at 9.7 fim has also been observed towards several galax¬ 
ies i ncluding LINERS, Seyfert galaxies, ULIRGs, and QSOs 
fe.g. IStu rm et ah l200fil: iHa o et al.ll 200Rl : ISpoon et al.l l2007l : 
ISmith et aI.ll201C : iMason et al.l l2012l l and can be used to 
constrain th e geometry and st ructure of the emitting nu¬ 
clear region (|Mason et al.ll2009l 'l. 

M31 with its proximity (785 ± 25 kpc; 

iMcGonnachie et al.l 1200^ 1 and rich observational databases 
provides the most detailed view of a star forming galaxy 
similar to the Milky Way. The acti ve star forming ring 
visible in 8 /rm Spitzer/IRAC images llBarmbv et al.l 120061 1 
provides evidence of abundant PAHs in M31. However, 
ISOC AM spectro-imaging observations dCesarskv et al.l 
Il998h of four regions in M31 including the nucleus and 
bulge of this galaxy showed very odd PAH spectra, bright 
at 11.2 fim with weak or no 6.2, 7.7, and 8.6 fim bands. 
Investigating this unusual PAH emission was the main 
motivation for the work described in this paper. 

The centre of M31 has a complicated physical struc¬ 
ture. It hosts a very inactiv e supermassive black hole with a 
mass of 0.7— 1.4 x 10® Mq llBacon et al.ll200ll : lBender et al.l 
l2005h and also has a lopsided nuc lear disk with two stel- 
lar components llLauer et al.l Il993l l and an A-star cluster 
llBender et al.l[2005l ~l. While M31’s nucleus is known to be 
inactive from an X-ray perspective llLi et al. I I 2 OIII I , mid- 
infrared indicators of its nuclear activity, such as infrared 



F„(8 nm)/F„(24 urn) 

Figure 1. 8 — 24/70 — 160 fira colour-colour diagram of M31 
obtained from IRAC and MIPS imaging. Colour-coding of points 
represents 24 /rm flux, from faintest (blue, dark grey) to average 
(green, light grey) to brightest (red, medium grey). The plot is 
divided into 9 regions (black grid), and the observations were 
made to cover those regions subject to a 24 fim brightness cut. 
The black triangles indicate colours of the regions we observed. 


excess or spectral features of silicates, have received rela¬ 
tively little attention. The higher spatial resolution available 
in observations of a very nearby galaxy like M31, compared 
to lumin ous, distant object s such as ultra-lumi nous infrared 
galax ies dSpoon et al.ll20()7h or nearby Seyferts llMason et al.l 
l2009l) . makes exploring its mid-infrared spectrum worth¬ 
while. 

We employed mid-infrared spectral maps from the 
Spitzer/lufTaied Spectrograph (IRS) from 12 regions of M31 
for a further investigation of its infrared properties. This 
sample includes the nucleus, bulge, an active region in the 
star-forming ring (all previously observed by ISOCAM), and 
9 other regions chosen to cover a range of properties as 
described in Section o We also obtained the highly pro¬ 
cessed version of IS OGAM observations of M31 and com¬ 
pare them with the IRS results in Section \M\ Section \n\ 
discusses PAH intensity ratios. In Section 14.2[ we investi¬ 
gate the relationship between PAH equivalent widths, radi- 
ati on hardness, and metalli city and compare to that found 
by lEngelbracht et all (l2008l ) and I Gordon et al.l (I 2 OO 8 I ). Sec- 
tionl^discusses the dust properties of the nucleus. The paper 
is summarized in Section |6l 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 IRS observations 

We obtained mid-infrared spectral ma ps of 12 regions in 
M31 using the Spitzer/IRS instrument (iHouck et al.l[2003) 
covering wavelengths from 5 to 21 microns. A background 
observation was also made off the galaxy along the minor 
axis for use in background subtraction from the data cubes. 
The 12 target regions include the nucleus, the ‘bulge’ and 
‘active’ regions previously observed by ISOGAM (the lat¬ 
ter is Region 9 in our sample), and 9 other regions selected 
to cover a range of metallicities and dust temperaturesQ 


^ One additional spectral map in M31 is available in the Spitzer 
archive (AOR key 12019200); unfortunately it does not cover the 
5-13 fim region, which is the major focus of this paper. We there¬ 
fore elected to not include these observations. 
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Mid-infrared spectroscopy of MSI 3 



Figure 2. An 8 micron negative IRAC image of M31 (dark colours indicate higher flux. iBarmbv et al.ll2006l ~). Small white rectangl es 
(30" X 50") show the regions that we observed, and larger squares (192" X 192") show the regions observed bv ICesarskv et ^ lll998l '). 


These 9 regions were chosen by c onvolving the IRAC 8 pm 
llBarmbv et al.l[200^ 1 and MIPS (iGordon et al.l[200^ maps 
to the same resolution and constructing an 8 — 24/70 — 160 
pm colour-colour diagram (Figure[T]). This colour space was 
used to give a rough definition of the types of spectral en¬ 
ergy distribution; the dense region in the plot was split into 
a 3x3 grid and one pixel in each grid region (subject to a 
24pm brightness cut) was selected for spectroscopy. The lo¬ 
cations of the observed regions are shown in Figure [51 and 
their coordinates and metallicities are given in Table [T] Ex¬ 
cept for regions 5 and 8, all of our mapped regions con¬ 
tain an H II regio n with an opti c al spe ctroscopic metallicity 
measurement by I Sanders et al.l (l2012l): Table [T gives the 
measurements f rom the method ISanders et ahl l 20121 ) de¬ 
note “NOG N2” llNagao et al.ll200^ ). For regions 5 and 8 we 
estimated metalliciti es using the M31 radial gradient that 
ISanders et alJ (l2012l) fit to the Hll region “NOG N2” mea¬ 
surements: 12 + log [0/H] = 9.09 — 0.00195i?gc. 

For our observations we used the IRS Short-Low (SL) 
and Long-Low (LL) modules, which cover wavelengths from 
5 to 21 microns. These modules have resolving power in the 
range 60-130. Each low-resolution module is divided into 
two sub-slits which provide spectroscopy in either first or 
second order. They are denoted as SLl (7.5-14.5 pm), SL2 
(5.2-7.6 pm), LLl (20.5-38.5 pm, not used in these obser¬ 
vations), and LL2 (14.5-20.75 pm). The data were obtained 
in ‘no peak-up’ mode resulting in a positional uncertainty 
of 1". All M31 regions were observed in September 2007 as 
part of G. G. Fazio’s Guaranteed Time (program ID 40032). 
The map size was based on the size of the IRS slits (SL: 
3.6" X 57", LL: 10.5" x 168"). Each region was covered by 
18 overlapping observations of the SL slit and 11 overlapping 
observations of the LL slit making the map size 32" x 57" for 
SL and 58" x 168" for LL. Figure [3] shows an example of the 
slit arrangement. For the brighter regions (nucleus, bulge), 
ramp times of 14 s (SL) and 30 s (LL) were used, while for 
the fainter regions, ramp times of 60 and 120 s were used 
respectively. Background observations were taken with each 
module (2 per ramp time). Because all of the targets are in 
the same part of the sky, a common background observation 
was used for multiple targets to subtract the background 
emission. 


Intensity (MJy/sr) 

10 20 30 40 so 60 70 00 90 



Figure 3. The 7.6 pm plane constructed from the SLl data cube 
of the nucleus, showing the arrangement of slits used to cover the 
region. Black boxes outline the footprint of the SLl maps (the off- 
centre SLl map is from observations made when SL2 was centred) 
and the green box outlines the LL2 map. Blue and red slits show 
how each map was covered using overlapping slit positions. 


2.2 IRS Data Reduction 

The data were reduced through the SSC pipeline (ver. 
S17.2.0), and the maps were assembled using the GUBISM 
program (ISmith et al.l l2007al ). Bad pixel removal was also 
done using GUBISM, and the background observations were 
used to subtract the background emi ssion from the s e cube s 
following the method outlined by [Gordon et al.l ll2008l) . 
Spectra were extracted using a 30" x 50" rectangular aper¬ 
ture, which corresponds to 114 x 190 pc at the distance of 
M31. The aperture size was selected to cover the overlap¬ 
ping area of the SL and LL modes; all the IRS maps cover 
more area than considered here. In order to study the spatial 
variation of the emission near the nucleus, we also extracted 
spectra from two smaller regions within that map; these will 
be further discussed in Section^ The spectrum of NGC 206 
is very noisy and was removed from our analysis. 

There is wavelength overlap between the SLl and SL2 
spectra and also between the SLl and LL2 spectra. To gen- 
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Table 1. Spitzer/IRS Target Locations in M31 


Name 

R.A. (J2000) 

Deck (J2000) 

Rgc*' 

12 -t log(0/H)° 

Nucleus^ 

00*'42“44?31 

41°16'09'.'4 

0.0 


Bulge“ 

00*'42™35!00 

41°21'0T.'0 

4.7 

8.90 ±0.03 

Region 1 

00*'41™30!41 

40°43'0778 

12.4 

9.20 ±0.20 

Region 2 

00*'45™22!85 

41°38'5371 

13.0 

9.07 ±0.02 

Region 3 

00''40™37;37 

41°01'2974 

12.1 

8.85 ±0.01 

Region 4 

00‘'41™17!86 

41°07'0978 

8.7 

8.89 ±0.06 

Region 5 

00‘'43™39f57 

41°19'0371 

7.0 

8.93 ±0.08 

Region 6 

00''43™35!72 

41° 23'1570 

4.3 

8.73 ±0.08 

Region 7 

00*'40™53!98 

40°58'5879 

8.7 

8.40 ±0.08 

Region 8 

00‘'42™21!60 

41°07'1774 

3.1 

8.94 ±0.08 

Region 9*^ 

00*'41™00;00 

40°36'20'.'3 

13.5 

8.86 ±0.02 

NGC 206 

00‘'40™20!20 

40°44'54'.'0 

9.8 


Background 

00‘'44™41!80 

40°58'5670 

29.5 



“Regions with ISOCAM data. 

^De-projected gala ctocentric distance in kpc. 

“Metallicities from [Sanders et al.l ll2012h . except for Regions 5 and 8 
where metallicities are estimated from the radial metallicity profile. 


erate a single spectrum for each M31 region it is necessary 
to combine the spectra and account for photometric offsets 
between them. Such offsets are commonly seen in IRS spec¬ 
tra extracted from extended regions. For most regions, the 
SLl and SL2 flux densities were quite well-matched over the 
wavelength overlap region (7.5 < A < 7.6/rm), and those 
two orders were joined without offset. In two cases with off¬ 
sets between SLl and SL2, SL2 and the ‘bonus order’ SL3 
were well-matched. We therefore combined the SLl and SL3 
spectra by adding a constant to the SLl spectra so that 
they matched the SL3 average. The SLl and SL2 orders 
were then combined as described above. After this proce¬ 
dure there was still a noticeable mis-match between the 
SL and LL spectra. We addressed th is by scaling the SL 
spectra to match IRAC 8 /rm fluxes dBarmbv et al.ll2006l l 
in the same apertures used to extract the IRS spectra. An 
extended source aperture correction of 0.824 was applied to 
the IRAC fluxes; this value wa s computed using the formula 
in the IRAC Data Handbook (ISoitzer Science Centerir2013l l 
with the radius of a circular aperture having the same area 
as the extract ion aperture. The Spitzer sy nthetic photom¬ 
etry software (ISoitzer Science Center! [20121 ') was then used 
to quantify the colour correction for each spectrum, i.e. the 
multiplicative factor K between the IRAC photometry over 
its broad bandpass and the IRS flux density at the cen¬ 
tre of the bandpass. For each spectrum we computed the 
additive offset between the colour-corrected IRS flux den¬ 
sity and the IRAC aperture-corrected flux density, and used 
thi s to correct the SL spe ctra; this method was also used 
by ISandstrom et al.l (|2012| L The SL spectra corrected with 
additive offsets were closely-matched to the LL spectra. Al¬ 
though a multiplicative offset could also have been used to 
match the IRS photometry to IRAC, the resulting change 
in spectral slope increased the size of the SL/LL offset. 


2.3 ISOCAM Data Reduction 

To compare our results with those of ICesarskv et ahl lll998ll . 
we retrieved the highly processed ISOCAM data provided by 
iBoulanger et al.l (l2005l l for the nucleus, bulge, and region 9. 


Intensity (MJy/sr) 

0 IS 30 45 60 75 90 105 

: :ia: 


$l 16'00' 



Figure 4. 11.2 fim negative image of the ISOCAM data cube 
from the nucleus of M31. The black box shows the size of the 
aperture (30'' x 50") used to extract IRS spectra. 


The ISOCAM data were obtained with circular variable fil¬ 
ters (CVFs) over a 3' x 3' field of view at a scale of 6" per 
pixel. The wavelength range covered was 5.15-16.5 /rm at 
a resolutio n of A/AA ss 45 ; the ISOCAM instrument is de¬ 
scribed bv ICesarskv et all lll996h . An example image of the 
ISOCAM data is shown in Figure [I] For the three regions, 
we extracted spectra using the same 30” x 50” aperture as 
for the IRS data. 


3 SPECTRAL CHARACTERISTICS AND 
FEATURE MEASUREMENTS 

The final processed IRS spectra are shown in Figure [S] ex¬ 
cept for the nucleus, discussed in Section [5] All of the main 
PAH features, including the 6.2, 7.7, 8.6 and 11.2 /rm bands, 
are clearly visible for all the regions shown here. The IRS 
spectra also show atomic line emission such as [Ar ll], [Ar ill], 
[S III], [S iv], [Ne ll], [Ne lll] and molecular H 2 emission at 
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Figure 5. Observed IRS spectra and detailed PAHFIT decompositions (see Section 13.211 . Black circles show the observed data, and 
green, blue, orange, magenta and red lines represent the modelled dust continua, PAH features, atomic and molecular lines, starlight 
continuum and the fit respectively. The grey line shows the total modelled continuum. Vertical scales differ in different panels. Spectra 
from the nucleus and NGC 206 are not shown here. 
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12.3 ^m. Some of the spectra display a contribution to the 
continuum from starlight emission. 


3.1 ISOCAM versus IRS 

As mentioned in the Introductio n, based on ISOCAM ob¬ 
servations of four regions in M31. [^sarskv et al.l (Il998ll re¬ 
ported suppression of the 6-8 /rm PAH bands and a strong 

11.2 /rm PAH band. The 11.2 /rm band profile was found to 
vary from region to region. In addition. IPagani et al.l (Il999l l 
confirmed that the star-forming ring in M31 seemed to show 
very weak PAH emission in the 6 to 8 fj,m region. However, 
the IRS spectra presented here in Figure [5] do not show 
such unusual behaviour. Indeed all regions, except the nu¬ 
cleus, show a normal mid-infrared spectrum similar to other 
nearby star-forming galaxies. 

Until 2005, ISOCAM-CVF data were not properly back¬ 
ground subtracted, and they were contaminated with zo¬ 
diacal emission and stray light. Differential spectra be¬ 
tween regions of relatively strong and weak emission were 
previously used to overcome this proble m (more details 
abou t the differential sp ectra are given by ICesarskv et al.l 
Il998l . IPagani et al.lll999h . In 2005, all ISOCAM-CVF data 
were reprocessed and corrected for the zodiacal emission 
llBoulanger et al.]l2005l l. We obtained these newly-processed 
ISOCAM spectra from three regions in our IRS sample (see 
Section [2.311 in order to compare them with the correspond¬ 
ing IRS spectra. Figure |6] shows that although the relative 
feature intensities in the IRS and ISOCAM spectra differ 
in detail, the spectral shapes are almost identical; the re¬ 
processed ISOCAM data do not agree with the previous dif¬ 
ferential spectra. Neither the bulge nor Region 9 shows any 
depletion in 6 to 8 /rm features. The new ISOCAM reduc¬ 
tion appears to eliminate the discrepancies between ISO¬ 
CAM and IRS, resulting in less ‘strange’ infrared spectra 
for M31. For the remainder of this work, we discuss only the 
IRS spectra. 

3.2 PAHFIT 

The PAH features in IRS spectra are often blended with 
neighbouring PAH features and atomic lines. Therefore mea¬ 
suring the strength of PAH features is difhcult . To achieve 
this tas k a tool called PAHFIT, introduced bv ISmith et al.l 
ll2007bll . was used. PAHFIT is an IDL based tool designed 
for decomposing Spitzer IRS low-resolution spectra of PAH 
emission sources. PAHFIT uses six main components to ht 
the surface brightness. These are starlight continuum, fea¬ 
tureless thermal dust continuum, pure rotational lines of H 2 , 
fine-structure lines, dust emission features and dust extinc¬ 
tion. The starlight is represented by blackbody emission at 
a fixed temperature of 5000 K, and the dust continuum is 
represented by 8 modified blackbodies (emissivity propor¬ 
tional to i/^) at fixed temperatures of 35, 40, 50, 65, 90, 
135, 200, and 300 K. The final fit obtained with PAHFIT 
does not necessarily use all eight dust components. The in¬ 
frared extinction is considered as a combination of a power 
law plus silicate absorption features with peaks at 9.7 and 
18 pm. Line features are represented by Gaussian profiles 
with widths set by the instrumental resolution and dust fea¬ 
tures are represented by Prude profiles; m ore details about 
PAHFIT are given bv ISmith et al.l ll2007bll . 



Wavelength (pm) 


Figure 6. Comparison of IRS and re-processed ISOCAM spectra 
for the Nucleus (top). Bulge (middle) and Region 9 (bottom) in 


M31. 


Initial attempts at htting the spectra with PAHFIT 
showed that some components were negligible, and to avoid 
over-fitting we re-ran the fits without these components. 
None of the IRS spectra shows significant silicate absorp¬ 
tion around 9.7 or 18 pm, and the extinction calculated by 
PAHFIT was almost zero for all the initial fits. Except for 
four regions (the bulge. Region 5, Region 6 and Region 8), 
the starlight contribution is also negligible. This is consistent 
with the fact that regions 5, 6, and 8 are the innermost of 
the regions. We adjusted the PAHFIT input parameters to 
fix extinction to zero for all regions and starlight to zero for 
all but the four regions above. Typically only two or three 
thermal dust components had significant power in our fits, 
but we did not fix the unused components to zero. Regions 3 
and 9 were found to have very low dust continuum emission 
compared to other spectra, possibly because of noisy data at 
short wavelengths. However the other features in these spec¬ 
tra appear to be fit correctly. The spectrum of the nucleus 
shows silicate emission, which is not included in PAHFIT; 
our modifications of PAHFIT to analyze this spectrum are 
discussed in Section [S] 


3.3 PAH features 


PAHFIT returns fluxes and equivalent widths (EQWs) of 
PAH features, which are given in Tables [5] and [21 For easier 
comparison with previous work, the tables give measure¬ 
ments for the 7.7, 11.2, 12 .7 and 17.0 pm PAH complexes 
as defined bv [Smith et al.l ll2007bf) . as computed from the 
individual constituent features measured by PAHFIT. The 
intensities of the features do not include any contribution 
from the continuum, but the equivalent width computed by 


EQW 


■dX, 


7i/,c 


( 1 ) 
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Table 2. PAH Emission Line Strengths" 


Region 

5.7 pbTo. 

6.2 pLm 

7.7 

b 

pLxn. 

8.3 /4m 

8.6 /4m 

10.7 , 

(4m 

11.: 

•» h 

1 /4m 

12.0 ; 

um 

12.7 um'’ 

17.0 

1 /4m^ 

Bulge 




1.32 

± 

0.08 

7.7 

± 

0.6 

1.1 

± 

0.2 

0.7 

± 

0.1 

0.07 

± 

0.03 

1.85 

± 

0.08 

0.49 

± 

0.05 

1.02 

± 

0.09 

1.39 

± 

0.04 

Region 1 

0.36 


0.05 

1.20 

± 

0.04 

3.8 

± 

0.2 

0.46 

± 

0.04 

0.50 

± 

0.03 

0.08 

± 

0.01 

1.18 

± 

0.02 

0.32 

± 

0.02 

0.54 

± 

0.02 

0.58 

± 

0.02 

Region 2 

0.27 


0.03 

1.10 

± 

0.03 

3.7 

± 

0.2 

0.33 

± 

0.04 

0.70 


0.03 

0.053 

± 

0.008 

1.13 

± 

0.02 

0.23 


0.01 

0.51 

± 

0.03 

0.51 

± 

0.02 

Region 3 

0.3 

± 

0.1 

0.9 

± 

0.1 

3.9 

± 

0.6 

0.7 

± 

0.1 

0.3 

± 

0.1 




0.68 

± 

0.07 

0.21 

± 

0.05 

0.47 

± 

0.07 

0.45 

± 

0.04 

Region 4 

0.14 


0.04 

0.56 

± 

0.03 

2.1 

± 

0.2 

0.26 

± 

0.04 

0.41 

± 

0.03 

0.029 


0.008 

0.70 

± 

0.02 

0.12 

± 

0.01 

0.31 


0.02 

0.44 


0.02 

Region 5 




0.24 

± 

0.04 

0.79 

± 

0.07 

0.12 

± 

0.04 

0.21 

± 

0.03 

0.032 

± 

0.008 

0.45 

± 

0.02 

0.09 

± 

0.01 

0.22 

± 

0.01 

0.33 

± 

0.04 

Region 6 




0.26 

± 

0.03 

0.8 

± 

0.2 

0.10 

± 

0.03 

0.13 

± 

0.03 

0.029 

± 

0.008 

0.38 

± 

0.02 

0.07 

± 

0.01 

0.16 

± 

0.01 

0.29 


0.02 

Region 7 

0.21 

zh 

0.03 

0.62 

± 

0.03 

2.0 

± 

0.2 

0.30 

± 

0.04 

0.45 

± 

0.03 

0.058 

± 

0.008 

0.77 

± 

0.02 

0.18 

± 

0.01 

0.37 


0.03 

0.47 


0.04 

Region 8 




0.09 

± 

0.04 

0.22 

± 

0.07 

0.09 

± 

0.04 

0.10 

± 

0.03 

0.051 

± 

0.009 

0.16 

± 

0.02 







0.15 


0.02 

Region 9 




1.3 

± 

0.1 

4.4 

± 

0.6 

0.9 

± 

0.1 

0.5 

± 

0.1 

0.09 

± 

0.03 

1.32 

± 

0.07 

0.51 

± 

0.05 

0.83 

± 

0.07 

0.6 


0.1 


“Units are 10 W m Uncertainties are as calculated by PAHFIT (see text for details). 
^PAH complexes, as described in text. 


Table 3. PAH Emission Line Equivalent Widths" 


Region 

5.7 /4m 

6.2 /4m 

7.7 Aim'’ 

8.3 /4m 

8.6 /4m 

10.7 /4m 

11.2 Aim'’ 

12.0 /4m 

12.7 Aim'’ 

17.0 Aim'’ 

Bulge 


1.2 ± 0.2 

4.1 ± 0.4 

0.51 ± 0.07 

0.30 ± 0.06 

0.03 ± 0.01 

0.78 ± 0.04 

0.22 ± 0.02 

0.49 ± 0.05 

1.16 ± 0.04 

Region 1 

0.39 ± 0.08 

1.2 ± 0.1 

3.4 ± 0.3 

0.43 ± 0.04 

0.47 ± 0.04 

0.09 ± 0.01 

1.45 ± 0.04 

0.43 ± 0.03 

0.76 ± 0.04 

1.26 ± 0.05 

Region 2 

0.28 ± 0.04 

1.02 ± 0.06 

3.4 ± 0.2 

0.32 ± 0.04 

0.70 ± 0.04 

0.07 ± 0.01 

1.58 ± 0.04 

0.35 ± 0.02 

0.85 ± 0.05 

1.33 ± 0.06 

Region 3“ 

4.3 ± 2.3 

8.3 ± 2.3 

18.5 ± 4.2 

2.8 ± 0.8 

0.9 ± 0.5 


2.1 ± 0.3 

0.7 ± 0.2 

1.6 ± 0.3 

2.4 ± 0.3 

Region 4 

0.28 ± 0.09 

1.0 ± 0.1 

3.7 ± 0.4 

0.48 ± 0.07 

0.77 ± 0.07 

0.07 ± 0.02 

1.67 ± 0.06 

0.31 ± 0.04 

0.82 ± 0.06 

1.65 ± 0.08 

Region 5 


0.12 ± 0.03 

0.61 ± 0.07 

0.10 ± 0.03 

0.20 ± 0.03 

0.05 ± 0.01 

0.77 ± 0.04 

0.17 ± 0.03 

0.50 ± 0.04 

1.6 ± 0.2 

Region 6 


0.10 ± 0.04 

0.6 ± 0.2 

0.10 ± 0.04 

0.14 ± 0.03 

0.05 ± 0.01 

0.77 ± 0.05 

0.15 ± 0.03 

0.42 ± 0.05 

1.7 ± 0.1 

Region 7 

0.32 ± 0.05 

0.86 ± 0.06 

2.2 ± 0.1 

0.44 ± 0.06 

0.69 ± 0.06 

0.12 ± 0.02 

1.81 ± 0.07 

0.48 ± 0.05 

1.08 ± 0.09 

2.2 ± 0.3 

Region 8 



0.10 ± 0.05 

0.09 ± 0.04 

0.10 ± 0.03 

0.09 ± 0.02 

0.30 ± 0.04 



0.62 ± 0.08 

Region 9“ 


(2.4 ± 1.1) X 10^ 

(1.4 ± 0.4) X lO'^ 

15.5 ± 5.6 

7.9 ± 2.8 

0.5 ± 0.2 

6.8 ± 1.1 

2.3 ± 0.6 

3.5 ± 0.8 

2.2 ± 0.8 


“Units are p.m. Uncertainties are calculated using a Monte-Carlo method (see text for details). 

^PAH complexes, as described in text. 

“Continuum for Regions 3 and 9 is very weak. Equivalent widths are highly uncertain and not considered in the analysis (see Section. 


is a measure of the ratio of the continuum emission (/i.,cont) 
to the line strength (/i.,feature = Iv — /j^.cont)- The contin¬ 
uum emission is mainly from dust grains much larger than 
PAH molecules. Hence, by studying EQWs of PAHs, we can 
study how the PAHs compete with the dust grains in the 
mid-infrared wavelengths. Uncertainties for line strengths 
are as calculated by PAHFIT’s propagation of statisti¬ 
cal uncertainties in the background-subtracted datacubes. 
Systematics and degeneracies in the fitting procedure are 
not included , and m ay dominate at higher signal-to-noise 
dSmith et al.ll2007lJl . Uncertainties for the EQWs for each 
PAH feature are calculated using a Monte-Carlo method. In 
that method, for each region, PAHFIT was run 500 times on 
randomly generated data points normally distributed within 
the uncertainties of the spectrum. PAHFIT returned 500 
EQWs for each PAH feature, and the standard deviation 
of EQWs for a given feature was taken as its uncertainty. 
The EQWs from Regions 3 and 9 were removed from fur¬ 
ther analysis because the negligible dust continuum for these 
spectra makes the EQWs highly uncertain. 



Figure 7. Log([Ne lll]/[Ne ll]) vs log([S lv]/[S III]) for the M31 
regions in our sample (b lack dots) and the starburst sample from 
lEngelbracht et al.l ll2008ll (open dots). The straight line is the line 
of best fit for the starburst sample. Arrows show 3 (t upper limits. 


3.4 Atomic line features 

PAHFIT also returns the line strengths and uncertainties for 
atomic lines, listed in Table 01 Not all lines were detected 
by PAHFIT, so we calculated upper limits for non-detected 
linesQ Line ratios of [Ne iii]/[Ne ii] and [S iv]/[S iii] 18 


hav e been used as an i ndica tion of the radiation hardness, 
and lEngelbracht et al.l (|2008h defined a combination of these 
two line ratios as a ‘radiation hardness index (RHI)’: 


RHI = 



[Ne III] 
[Ne ii] 


-b 


0.71 -b 1.58 log 


[S IV] 

[S III] 


/2 


( 2 ) 


given by PAHFIT. The peak intensity was taken to be 3 times 
2 To find the upper limits for the flux of missing atomic lines, the RMS, where RMS is the root mean square of the noise at the 

we assumed the line to be a Gaussian profile with a FWHM as position of a missing line. 
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Table 4. Atomic Emission Line Strengths“ 


Region 

[Ar ii] 

7.0 Aim 

[Ar III] 

9.0 /xm 

[S IV] 

10.5 Aim 

[Ne ll] 

12.8 fim 

[Ne III] 

15.5 fim 

[S III] 

18.7 /xiri 

Bulge 

< 0.12 

0.17 ±0.03 

< 0.11 

0.07 ±0.01 

0.025 ± 0.006 

0.007 ± 0.003 

Region 1 

< 0.05 

< 0.06 

0.020 ± 0.004 

0.020 ± 0.004 

< 0.01 

0.008 ± 0.001 

Region 2 

< 0.05 

< 0.06 

< 0.02 

0.022 ± 0.004 

< 0.01 

0.003 ±0.002 

Region 3 

< 0.15 

0.09 ± 0.02 

< 0.10 

0.03 ±0.01 

0.020 ± 0.004 

0.015 ±0.003 

Region 4 

< 0.04 

< 0.04 

< 0.02 

< 0.02 

< 0.01 

0.004 ±0.002 

Region 5 

< 0.04 

< 0.02 

< 0.02 

< 0.01 

0.009 ± 0.004 

0.008 ±0.004 

Region 6 

0.02 ±0.01 

0.014 ±0.007 

< 0.01 

< 0.01 

0.019 ±0.002 

0.019 ±0.002 

Region 7 

< 0.04 

< 0.04 

0.008 ± 0.003 

0.035 ± 0.005 

< 0.01 

0.027 ± 0.004 

Region 8 

< 0.04 

0.017 ±0.006 

< 0.02 

< 0.01 

0.041 ± 0.002 

0.023 ±0.002 

Region 9 

0.09 ± 0.04 

0.12 ±0.03 

< 0.09 

0.13 ±0.01 

< 0.04 

0.05 ±0.02 


“Units are 10 W m Uncertainties are as calculated by PAHFIT (see 
text for details). Upper limits (3cr) are indicated with a < mark. 


Here, 1.58 and 0.71 are the slope and intercept of the 
[Ne lll]/[Ne ll] vs [S lv1/[S lll] plot (Figure [71) for the star- 
burst sample from lEngelbracht et al.l (l2008l) . Figure [7] com¬ 
pares the atomic line emission from the selected regions of 
M31 to the starburst galaxy sample; although none of our 
spectra has detections of all four lines, our limits are mostly 
consistent with the trend. To compute RHI in the case of 
missing lines, we used the term in Equation[2]for which both 
lines are detected. For regions with only one line detected per 
element, we computed RHI upper limits using Equation [2] 


4 RESULTS: EXTRANUCLEAR REGIONS 
4.1 PAH band ratios 


Both the 6.2 and 7.7 fim features are thought to be coming 
from ionized PAHs and the 11.2 fim feature from neutral 
PAHs. Therefore we expect to see a correlation between the 
intensities of 6.2 and 7.7 /rm PAH features normalized by the 
11.2 fim feature. Figure [S] compares the PAH flux ratios of 
7.7/11.2 and 6.2/11.2 features. The figure shows a good cor¬ 
relation between these two PAH line ratios, consistent wit h 
that of the SINGS sample shown bv ISmit h et al.| (|2007bll. 
A sim ilar correlation was also reported by Galliano et al.l 
(l2008l) for a sample of galaxies and a ha ndful of extended 
Hll regions and by IVermeii et al.l (l2002l ) for Galactic and 
Magellanic Gloud Hll regions. This provides evidence that 
the PAH emission from M31 is not unusual. 


4.2 PAH equivalent widths versus radiation 
hardness and metallicity 

As mentioned in the introduction, PAH equivalent widths 
tend to decrease as radiation hardness increases, and as 
metallicity decreases llCalzetti et al.l[2010ll . 

Figure [9] shows the equivalent widths of the 8 ^m fea¬ 
ture (a combination of the 7.7, 8.3 and 8.6 ^m PAHFIT 
components) and the 7.7 and 11.2 ^m features as a func¬ 
tion of both metallicity and RHI . The starburst galaxy sam- 
ple of lE ngelbra cht et al] ll2008l) and Hll regions in MlOl 
dGordon et al.ll2008l l are plotted for comparison. The M31 
regions cover approximately the same range of measure¬ 
ments as the starburst galaxies and MlOl Hll regions; the 



Log(PAH_6.2/PAH_11.2) 


Figure 8. Ratios of PAH feature fluxes (7.7 /rm/11.2 fim versus 
6.2 ^m/11.2 fim) for 10 regions in M31. Open squares represent 
the centra l regions of nearby g alaxies as observed in the SINGS 
sample bv ISmith et al.l ll2007bh . 


outliers are regions 3 and 9 with high equivalent widths (un¬ 
certain due to low continuum, see Section [3.211 and region 8 
with low equivalent widths (noisy spectrum as well as sub¬ 
stantial modelled contribution from starlight, see Figure [5]). 
No clear trend is defined by the M31 data, unsurprising given 
the uncertainties and the limited number of regions. We do 
not have enough data from low-metallicity regions in M31 to 
observe the expected decrease of PAH EQW with decreasing 
metallicity; however the M31 data occupy the same region of 
parameter space as the MlOl data and we conclude that the 
EQWs of the regions in M31 are consistent with previously 
published “normal” values. 


5 RESULTS: THE NUCLEAR REGION 

Examining the Spitzer-IRS spectral data cube for the nu¬ 
clear region, we noticed that different spectral features vary 
spatially within this region. The 11.2 pm PAH emission is 
discrete and patchy (Figure [RH top). Indeed, the majority 
of the 11.2 pm PAH emission is from a region centered at 
(00:42:43.947, -1-41:16:22.92), 15” north of the nucleus and 
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Figure 9. Equivalent width of PAH feature versus metallicity (left panels) and radiation hardness index (RHI; right panels). Equivalent 
widths are not normalized ; metallicities of th e M31 regions have had 0.35 dex subtracted to account for the offset between direct and 
strong-line measurements d Crolev et al.l l20l^ . Solid points: M31 sample; triangles represe nt upper limits . Open points: starburst galaxy 
sample from lEngelbracht et ahl (l2008l . bottom panels) or H II regions in MIDI observed by [Gordon et all 1 I 2 OO 8 I . middle and top panels). 


not from the nucleus itself. Weaker 11.2 /rm PAH emission is 
also found near the edge of the map peaking at 00:42:45.497, 
-h41:15:43.97 and near 00:42:45.869; -h41:16:11.38. The lo¬ 
cations of the two weaker 11.2 pm PAH emission peaks are 
also near positions exhibiting CO(2 -l) line emission (#36 
and 28 of iMelchior fc CombedBoiSl . the strongest 11.2 pm 
PAH emission peak is outside the CO FOV). On the other 
hand, the centre shows no PAH emission, but it does have 
clear silicate emission around 9.7 pm (Figure [TOl bottom) 0 
Although [Nelll] 15.5 pm line emission is strong at the three 
locations with 11.2 pm PAH emission and weak(er) at the 
nucleus, it is also present across the nuclear region and thus 
spatially distinct from the silicate and PAH emission. 

Some structures seen in the optical near the M31 nu¬ 
cleus, such as the double nucleus and cluster of young stars, 
cannot be resolved with the data discussed here. The IRS 
spectral cubes have pixel sizes of l”8, and the SL PSF 
FWHM is 2.5-3” depending on the wavelength, while the 


® The spatial resolution and pixel scale of the ISOCAM data are 
not sufficient to resolve the centre and north regions. 


two components of the nucleus are separated by only 0'.'5 
llBender et al.l[2005l) . However, the north 11.2 pm PAH emis¬ 
sion appear to be marginally spatially resolved: its radial 
profile has a full width at half maximum (FWHM) of 5-7” 
(corresponding to 19-27 pc). 

To examine this spatial variation towards the nucleus 
in more detail, we extracted spectra from the centre 
and the North regions using 9” x 9” square apertures 
in addition to the full 30” x 50” M31 nuclear spectrum 
(further referred to as nuclear spectrum, north spectrum 
and full spectrum respectively), as shown in Figures 1111 
and m Both spectra show a blue continuum and atomic 
fine-structure lines of Ne and S, but they exhibit distinct 
dust emission consistent with the spatial maps: PAH 
emission is detected in the north spectrum while silicate 
emission is dominating in the nuclear spectrum. In addition, 
weak H 2 emission is present at 17 pm in the north spectrum. 

To investigate the stellar contribution to these spectra, 
we made a comparison wit h the IR spectra of NU Pav, an 
M giant (ISloan et ahl 12013 ) and the templat e spectrum of 
quiescent elliptical galaxies as constructed bv iKaneda et ah] 
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RA U2000) 


Figure 10. Integrated strength of the 11.2 /j.m PAH emission (left) and the silicate emission (from 9 to 11 /im, continuum subtracted; 
right) around the nucleus of M31. The black boxes in the left panel are the 9" x 9^^ sub-apertures (centre and north region) used to extract 
spectra. The small white cir cle in the right pan el shows a position halfway between the two components of the nucleus (coordinates for 
PI and P2/P3 in Table 1 of iGarcia et aLll201(t) . These components are separated by 0^^5, which is well below the IRS spatial resolution. 
The contours in the right panel represent the distribution of 3.6 fim luminosity, indicating the orientation of the M31 bulge. 



Wavelength (//m) 

Figure 11. Mid-infrared spectrum of the nucleus of M31 (thick line; full spect rum) over-plotted with spectra extracted close to the nuclei 
of 6 nearby galaxies which have similar spectral shapes dSmith et al.|[^07blb The comparison spectra are offset by arbitrary amounts 
for visibility; they follow the same vertical order as in the legend. NGC 4552, NGC 1404 and NGC 4125 are elliptical galaxies and NGC 
4594 and NGC 2841 are spiral galaxies. NGC 1316 is a lenticular galaxy. The inset shows the spectra extracted from the centre region 
of the M31 nucleus and from the north region shown in Figure fTOl 
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Wavelength [/i.m] 


Figure 12. The spectra of the M31 nucleus (nuclear, north and 
full) compa red with the M giant N U Pav and the quiescent galaxy 
template of iKaneda et all l2008h . normalized to the 5.4 pm flux 
density. The legend gives the objects from lowest to highest flux 
at the longest wavelengths. 


ll2008l) to investigate the PAH emission in elliptical galaxies 
(in particular the 7.7 pm PAH band). The IR spectrum 
of NU Pav is very similar to this template spectrum, but 
it lacks silicate emission and dust continuum emission at 
the longer wavelengths. Assuming that the NIR fluxes are 
completely due to the stellar component, we can estimate 
the contribution of the stellar component to the 5.4 pm flux 
by normalizing the spectrum of NU Pav and the full M3 1 
spectrum to the H band flux (followingJVe£a^_et_^ |201^. 
NIR fluxes for NU Pav are taken from lGezari et aT 1 19991 ) 
and for M31 are determined by integrated photom etry of 
the 2MASS Large Galaxy Atlas Jjarrrtt et al'l[200^1 images 
over the spectral extraction regions. This allows us to deter¬ 
mine the stellar contribution to the 5.4 pm flux of the full 
M31 spectrum, which is found to be 96%. Therefore, in order 
to compare with the quiescent galaxy template, we further 
assume that 100% of the flux at 5.4 pm is stellar; Figure [T^ 
shows a comparison of the spectra when normalized at 5.4 
pm. The match between the spectra up to 8 pm is very 
close, with no evidence for non-stellar emission in the M31 
spectra. Furthermore, it is striking that the nuclear, full 
and north spectra of M31 match each other extremely well 
up to ~14 pm, including the 9 pm silicate emission, except 
for the PAH emission and fine-structure line emission. In 
contrast, the 9 pm silicate emission is less intense in the 
quiescent galaxy template and absent in NU Pav; both of 
these sources also lack PAH emission and Hne-structure 
line emission. Given that the presence and strength of 
the silicate emissio n, if stellar, depends on the age of the 
stellar population (IVillaume et al.ll2015l ). this resemblance 
suggests that the silicate emission is also stellar in nature. 
These authors predict the presence of silicate emission for 
stellar population ages of < 10®'® yr. M31 is known to have 
stars of this age group: a 200 Myr-old starburst in the 
central 1” plus an interme diate-age population both at the 
centr e and the inner bul ge jBender et al.l[2005l : ISaglia et al.l 
I 2 OIOI: iDong et al.l lioli) . Hence, the nuclear spectrum is 
largely dominated by the stellar component with little 
contribution from the nuclear region itself, which is mostly 
manifested by the fine-structure line emission. Indeed, while 


a weak dust continuum flux may be present at the longer 
wavelengths, the flux increase around 15 pm is likely 
due to the presence of the 18 pm silicate emission band. 
Therefore, assuming that the spectral shape of the stellar 
contribution does not change spatially within the M31 
nuclear region, we can use the nuclear spectrum as a tem¬ 
plate of the stellar contribution to the north and full spectra. 


The north spectrum exhibits clear PAH emission, in 
particular at 11.2 and 15-20 pm (Figure [12)). However, the 
typical 6-9 pm PAH emission is very weak or even absent. 
Such unusual PAH emi ssion has also been reported towards 
ellipt i cal galaxies (e . g. [X aneda et al.ll2 005l: iBregman et al.l 


I 2 OO 8 I : iKaneda et al.l I 2 OO 8 I : Vega et alJ 120101 ) and i 
low-lu m inosity active ga l actic nuclei (LLAGN, ISmith et al.l 
l2007bl ). iBregman et al.l (l2008l ) reported however that the 
shape of the stellar emission greatly influences the appear¬ 
ance, and thus apparent strength, of the 6-9 pm PAH emis¬ 
sion. Therefore, we apply two different corrections for the 
stellar component to the north spectrum in order to inves- 
tigate the ‘true’ PA H emission. First, we apply PAHFIT 
dSmith et al.ll200'7bh to the spectrum of the north region and 
rely on the stellar component within PAHFIT (a blackbody 
of 5000K) to account for the stellar contribution to the spec¬ 
trum (Figure fT^ left). This stellar component clearly does 
not take into account possible stellar spectral features such 
as silicate emission, H 2 O absorption (at 6.6 pm), SiO ab¬ 
sorption (at 8 pm) etc. Secondly, we first subtract the stellar 
contribution by using the nuclear spectrum (normalized at 
5.4 pm) and then apply PAHFIT to this stellar-subtracted 
north spectrum (Figure (TS) right). 

For both methods, the continuum emission is overesti¬ 
mated between ~ 8 to 11 pm indicating that the silicate 
emission may be slightly overestimated when using the nu¬ 
clear spectrum and confirming that a 5000 K blackbody (the 
stellar component in PAHFIT) is only a first-order approxi¬ 
mation for the stellar contribution. In both cases, we further 
subtract the total continuum emission as found by PAH¬ 
FIT and compare the residuals to a typical PAH emission 
spectrum as exemp lified by the Hll-type galaxy NGC 0337 
dSmith et al.ll200'7a ). Relative to the 11.2 pm PAH emission, 
the 15-20 pm emission is similar for both the north region 
and a typical PAH emission spectrum while the 6-9 pm PAH 
emission remains distinct in the north region. However, the 
6-9 pm PAH emission differs for the two methods applied. 
When using the stellar component of PAHFIT, a 6.2 pm 
PAH emission is present while no significant 7.7 and 8.6 pm 
emission is present. In contrast, when using the normalized 
nuclear spectrum as the stellar component, the 6.2 and 8.6 
pm PAH emission is largely absent, and the 7.7 pm PAH 
emission is weak. This is confirmed by the PAH intensity ra¬ 
tios found by PAHFIT: the 7.7/11.2 pm PAH intensity ratio 
is 0.79 for the M31 north spectrum whi ch is much lower tha n 
the average in the SINGS sample (3.6. ISmith et al1l2007bh . 
As mentioned, the two methods do not agree regarding the 
6.2 pm band, resulting in 6.2/11.2 pm PAH intensity ra¬ 
tios of 1.18 and 0.14 (for the PAHFIT and nuclear stellar 
component respecti vely) whereas the SINGS average is 1.1 
dSmith et al.ll2007bl ). In contrast, the 17/11.2 pm PAH in¬ 
tensity ratios are 0.46 and 0.52 (for the PAHFIT and nuclear 
stellar component respe ctively), which is co nsistent with the 
SINGS average of 0.53 JSmith et al.ll20079 ). 
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Figure 13. (top) PAHFIT result for the M31 north spectrum (left) and the M31 north spectrum subtracted by the stellar component 
(as defined by the M31 nuclear spectrum, right): fit (orange), continuum (magenta), individual dust components (blue), individual fine- 
structure lines and H 2 lines (green), (bottom) Continuum subtracted spectra of t he north region (left ) and the north region without the 
stellar component (right) compared with that of the Hll-type galaxy NGC0337 llSmith et al.ll^07bh . normalized to the peak intensity 
of the 11.2 /rm PAH emission. 


The nuclear component has the largest influence on 
the 7-9 /rm PAH emission due to the presence of silicate 
emission and SiO absorption while the 6.2 ^m PAH emission 
is less affected (due to weaker H 2 O absorption). Hence, 
it is remarkable that we recover a ’normal’ 6.2 ^m PAH 
emission by using a stellar blackbody while it is largely 
absent when a more ‘appropriate’ nuclear correction is 
applied. This testifies to the sensitivity of the deduced 
‘true’ PAH emission on the applied correction for the 
stellar component. Although the 7.7 ^m PAH emission is 
partly recovered after the stellar component correction, the 
7.7/11.2 PAH intensity ratio is still well below the average 
for the SINGS sample for both applied methods. Hence, 
until more accurate stellar templates are available, it is 
fair to conclude that the north spectrum does not show 
typical PAH emission . This is in contrast with the results of 
iBregman et al.l ll2008l) who recover a normal PAH spectrum 
for the elliptical galaxy NGC4697 after subtraction of a 
quie scent elliptical galax y spectrum. But it is consistent 
with iKaneda et ^ (l2008l ) who found that 80% of elliptical 
galaxies in their sample (of 15) remain anomalously weak 
7.7 ^m emission after correction for the stellar component. 

As mentioned earlier, similar atypical PA H emission has 
been detected towards elliptical galaxies (e.g. iKaneda et al.l 
I2OO8I : IVega et ^ l 2010 f ) and in some low-luminosi t y activ e 
galactic nuclei (LLAGN) reported bv ISmith et al.l l2007bl) . 
To place M31 in the context of these other galaxies, Figure 
compares the full 30” x 50” M31 nuclear spectrum with 
the nucle ar spectra of similar spectral shape from six SINGS 
galaxies (ISmith et al.ll2007bl )n All 6 SINGS galaxies share 


similar PAH feature characteristics to the M31 spectrum 
(e.g. weak 6-9 fim features), but none of them contains obvi¬ 
ous silicate emission. The SINGS galaxies with similar spec¬ 
tral shapes include three elliptical galaxies, two spirals, and a 
lenticular; there is some disagreement over the exact nuclear 
spectral types of these six galaxies llKennicutt et al.l l2003l : 


ISmith et al.ll20079 : iMoustakas et al.ll2010 ). All are classified 

as some form of LLAGN such as Seyfert or LINER (luminous 
AGNs were intentiona lly omitted from the SINGS sample; 
iKennicutt et al.ll2003l ). although they are by n o mean s the 
only LLAGNs in the SINGS sample. iLi et al.l ll2009t) con¬ 
cluded that the central black hole in M31 (MSI’*’) is currently 
inactiye, with direct obseryational signatures seen only at ra¬ 
dio and X-ray wayelengths, so finding additional signatures 
in the mid-infrared is of great interest. To our knowledge, no 
such signatures have been reported; bro adband mid-infrared 
imaging of the cent ral regions of M31 llDavidge et al.l[2006l : 
iBarmbv et al.ll2006l) did not identify unambiguous nuclear 
emission. The bluest part of the spectrum in Figure [TT] is 
dominated by the continuum, in agreement with the expec¬ 
tation that stellar light dominates at these wavelengths. 

Gould radiation from M31’*’ be responsible for the sup¬ 
pression of the 6-8 fim PAH f eatures compared to the 
11.3 nm fea t ure? A s discussed bv ISmith et al.l (l2007bl) and 
ISmith et al.l ll2010l) . inferring such a suppression must be 
done with caution, because the 6-8 ^m features are more 
susceptible to dilution by the stellar continuum. Several 
connections between PAH suppression and the presence of 
an AGN are possible, including destruction of small PAH 
molecules by a hard radiation field, modification of the struc- 


areas ranging from 2 to 8 kpc*^, whereas the M31 spectrum covers 
^ The IRS spectra for the SINGS galaxies were extracted over 0.02 kpc^. 
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ture of the PAH molecules, or weak ultraviolet continuum 
from lo w star formation rates leadiuR to decreased PAH exc i- 
tation (ISmith et aLll2007bl : iDiamond-Stanic fc Riekell201^ . 
In the latter case, the AGN is not the cause of the suppressed 
6-8 jj,m features but rather is only detectable because the 
nuclear star formation rate is low. Previous work has found 
low rates of star formation in the centre of M31: although 
iMelchior fc Combed ll2013h found a significant amount of 
cold gas in the centre of the galaxy, this gas does not ap- 
pear to be as sociated with current star formation (see also 
iLi et al.ll2009ll. In modelling the far-infrared spectral energy 
distribution, I Groves et al.l (l2012l l found that the old stellar 
population in the M31 bulge is sufficient to heat the observed 
dust; no young stellar population is needed. As shown, the 
nuclear spectrum is dominated by the stellar component and 
the dust continuum may be only weakly present at the longer 
wavelengths. This is consistent with a low star formation 
rate. We conclude that PAH feature ratios cannot provide 
direct evidence for radiation from M31*. 


6 SUMMARY AND CONCLUSIONS 

Spitzer/IRS spectral maps of 12 regions within M31 cover 
wavelengths 5-21 pm. The spectra from those regions, ex¬ 
cept for the nucleus, are similar to spectra obtained from 
other nearby star-forming galaxies. Early ISOCAM obser¬ 
vations towards 4 regions of M31 showing a suppression of 
the 6-8 pm features’ strength and a n enhancement of the 
11.3 pm feature intensity and FWHM llCesarskv et al.|[l998l l 
were likely affected by the incorrect background subtraction 
methods applied. Indeed, both the re-processed ISOCAM 
data and the IRS data (except for the nucleus) exhibit typ¬ 
ical PAH emission. 

We recover the well-known strong correlation between 
the 6.2/11.2 and 7.7/11.2 pm PAH intensity ratios that is 
governed by the PAH charge balance. This is another indi¬ 
cator that the PAH emission in all regions but the nucleus 
is typical. The PAH EQWs in M31 regions do not show a 
clear decreasing trend with increasing radiation hardness, 
but are consistent with previous results from other nearby 
galaxies. The distribution of PAH EQWs with metallicity 
is well within the r ange of the starburst galaxy sample of 
lEngelbracht et ^ (l2008l l. We did not have enough data 
from low-metallicity regions of M31 to observe the decreas¬ 
ing trend of EQWs at low metallicities which is visible in 
other galaxies. 

Different spectral features (11.2 pm PAH emission, 
[Nelll] 15.5 pm line emission, silicate emission) show distinct 
spatial distributions in the nuclear region. The mid-infrared 
spectrum from the nucleus of M31 shows a strong blue con¬ 
tinuum and clear silicate emission. This nuclear spectrum 
is largely dominated by the stellar component with little 
contribution from the nucleus itself. In addition to a blue 
continuum and silicate emission, the mid-infrared spectrum 
of a region north of the nucleus (15" off-nucleus) exhibits 
suppressed 6-8 pm PAH features and normal 11.3 and 17.0 
pm PAH emission. This atypical PAH emission is only seen 
elsewhere towards some elliptical galaxies and low lumi¬ 
nosity AGNs (e .g. ISmith et al.l l2007bl : iKaneda et al.l l2008l : 
IVeea et'aDl2010h . 
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